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Recent observations of the cosmic microwave background (CMB) anisotropies and the distribution of 
galaxies, galaxy clusters, and the Lyman a forest have constrained the shape of the power spectrum 
of matter fluctuations on large scales k < few h/Mpc. We explore a new technique to constrain the 
matter power spectrum on smaller scales, assuming the dark matter is a Weakly Interacting Massive 
Particle (WIMP) that annihilates at early epochs. Energy released by dark matter annihilation can 
modify the spectrum of CMB temperature fluctuations and thus CMB experiments such as Planck 
have been able to constrain the quantity /{(Jav) /m^ <1/88 picobarnxc/GeV, where / is the fraction 
of energy absorbed by gas, {aa,v} is the annihilation rate assumed constant, and is the particle 
mass. We assume the standard scale-invariant primordial matter power spectrum of Pprim(fc) ~ k^” 
at large scales k < kp, while we adopt the modified power law of Pprim(fc) ~ kp’‘{k/kp)'^^ at small 
scales. We then aim at deriving constraints on m^. For m-s > n^, the excess small-scale power results 
in a much larger number of nonlinear small mass halos, particularly at high redshifts. Dark matter 
annihilation in these halos releases sufficient energy to partially ionize the gas, and consequently 
modify the spectrum of CMB fluctuations. We show that the recent Planck data can already be 
used to constrain the power spectrum on small scales. For a simple model with an NFW profile 
with halo concentration parameter C200 = 5 and f {a^v) = 1/100 picobarnxc/GeV, we can limit 

the mass variance (Tmax 100 at the 95% confidence level, corresponding to a power law index 
mg < 1.43(1.63) for kp = 100 (1000) h/Mpc. Our results are also relevant to theories that feature 
a running spectral index. 

PACS numbers: 98.80.-k, 95.30.Sf, 98.62.Sb, 95.85.Ry 


I. INTRODUCTION 

The nature of dark matter is unknown and remains 
one of the greatest mysteries in astrophysics and cosmol¬ 
ogy. Weakly Interacting Massive Particles (WIMPs) are 
one of the leading candidates for the dark matter of the 
Universe, and large experiments are being conducted to 
detect WIMP dark matter through direct, indirect, and 
collider experiments. Some direct detection experiments 
such as DAMA [T] have observed an annual modulation 
consistent with the presence of dark matter particles of 
mass 8—15 GeV interacting with a spin-independent 
cross section of 0.01 — 0.1 femtobarn. Recent observations 
of the Milky Way center [2] by the Fermi gamma ray tele¬ 
scope also seem to indicate an excess of gamma rays, con¬ 
sistent with dark matter particles of mass m-^ = 31 — 40 
GeV annihilating at a rate (cTau) = (1.4 — 2.0) x 10“^® 
cm^/s. These exciting results are, however, inconsistent 
with the XENON [3] and LUX experiments [4], and are 
also disfavored by the non-detection of dark matter an¬ 
nihilation in the local dwarf galaxies 


‘Electronic address: anat01@me.com 


The cosmic microwave background has been shown to 
be an excellent probe of WIMP dark matter annihila¬ 
tion at high redshifts [iHiI]. Particle annihilation re¬ 
leases and injects energy into the cosmic diffuse gas, re¬ 
sulting in both ionization and heating. Free electrons 
scatter CMB photons and cause damping in the tem¬ 
perature anisotropy power spectrum at intermediate and 
small angular scales. Also the CMB polarization power 
spectrum is boosted at very large scales. Precise mea¬ 
surement of the CMB by Planck, WMAP, ACT, and 
SPT have already placed tight constraints on dark mat¬ 
ter properties [T^HTB] . Interestingly, the recent results 
from the Planck collaboration m constrain the anni¬ 
hilation parameter pann = /(^au)/??!^ < 3.4 x 10“^® 
cm®s“^GeV“^ = (1/88.3) pbxc/GeV, where / is the 
fraction of energy absorbed by gas. For realistic values 
of / « 0.35 for the bb channel |TH], and (aav) = 0.727 
pb X c [Hj , one obtains a bound on the dark matter mass 
> 22.4 GeV at the 95% confidence level[ii]. 

The CMB constraints on dark matter particle proper¬ 
ties are derived on the assumption that only annihilation 
of free ‘unbound’ particles contribute to the net energy 
release, i.e. annihilation of particles bound in nonlinear 
objects, “dark halos”, is not considered. The assump¬ 
tion is appropriate because, in the standard cosmology. 
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nonlinear objects appear only relatively late, at redshifts 
z ~ 30 — 50. Then the gas density is already small, and 
CMB photons do not interact with free electrons unless 
the gas is significantly ionized. Nonlinear halos would 
be, however, important if they formed much earlier, i.e. 
at redshifts z > 100. Such a case is possible if the pri¬ 
mordial density fluctuations have some excess power at 
small length scales. 

The simplest inflationary theories predict a nearly scale 
invariant primordial curvature power spectrum P-ji ~ 
resulting in a present day matter power spectrum: 

PUz,k) cx 

= ( 1 ) 

where T’prim(fc) cx fc”® is the primordial matter power 
spectrum on very large scales k < 10“^ /i/Mpc, T{k) 
is the transfer function, and D{z) is the growth factor. 
Data from the Wilkinson Microwave Anisotropy Probe 
(WMAP) and the Atacama Cosmology Telescope (ACT) 
was used by [20] to reconstruct the matter power spec¬ 
trum at wavenumbers 0.001 < k < 0.19 Mpc On 
smaller length scales, one may use the Sloan Digital 
Sky Survey measurements of the clustering of galaxies 
to probe scales up to fc ~ 0.2 h/Mpc. On even smaller 
scales, the flux power spectrum of the Lyman a forest 
may be used to probe the matter power spectrum for 
k < 2 Mpc“^. Unfortunately, at fc ^ 10 /i/Mpc, there 
exist no direct observations of the matter power spec¬ 
trum. Possible probes of the small scale power spectrum 
include the use of Type la supernova leasing dispersion 
HD, ultra compact mini halos P^25j . and the dissipa¬ 
tion of acoustic waves by Silk damping [26l [27] . 

In the present paper, we explore a new probe of the pri¬ 
mordial density fluctuations on very small scales. Let us 
consider a simple power-law for the matter power spec¬ 
trum at k > kp'. 


Uprini — Ak ^ k kp 

= Ak;^ {k/kpr^ k>kp, (2) 

which is consistent with all available observations pro¬ 
vided the pivot wavenumber kp is large enough, say 
kp > 10 /i/Mpc. Essentially, we examine if there is excess 
power on the relevant small length scales through energy 
injection from dark matter annihilation in the early uni¬ 
verse. The basic idea is as follows. The mass variance 
a^{z,M) is computed by integrating the dimensionless 
power spectrum A^(fc) = k^P(k)/2'K^ over a window 
function: 


tT^(z, M) 


D\z) 

£>2(0) 


dfc fc3p(fc) 

k 27r2 


W^{kR). 


( 3 ) 


the Planck cosmology, we find D{z = 0) = 0.757. The 
window function W(kR) for comoving scale R is conve¬ 
niently given by 


. 3 sin a; — a: cos xl 

Wix) = — - 5 - - 


( 4 ) 


and R is the comoving radius that encloses a mass M. 
Fig. [^a) shows the matter power spectrum for the stan¬ 
dard cosmology (solid lines, black) computed using the 
Eisenstein-Hu transfer function |28l [29]. Also plotted 
are curves for Ug = 0.96, and the modified power law of 
Eq. for rus = 1.1, 1.2, and 1.5. An exponential cutoff 
is imposed at the free streaming scale kfs ^ 10® h/Mpc 
[5nil32| . to account for the finite velocity dispersion of 
WIMP dark matter, and to make the integral in Eq. 
finite. This ensures that there is a minimum halo mass 
Almin ~ 1O“®M0. Panel (b) shows the standard devi¬ 
ation cr of density fluctuations, for these models. Note 
that CTniax = cr(Mmin) is very sensitive to mg, although 
as — 0.8 for all models. 


II. DARK MATTER ANNIHILATION IN HALOS 


Consider an overdensity of weakly interacting dark 
matter particles. The number of WIMPs in a volume 
SV is (p^/m^) SV, where is the density of WIMPs, 
and the probability of WIMP annihilation in a time St 
is {aav)St. The number of WIMP annihilations per unit 
time per unit volume is then equal to {a^v)p^^j Since 
each annihilation releases rriy^ of energy per particle, the 
total energy released per unit time per unit volume equals 


d£ _ (cTaU) 2 

dtdV rriyf. 


( 5 ) 


The energy per unit time due to particle annihilation 
in a bound halo of radius r 2 oo is obtained by integrating 
Eq.0 over the halo volume: 


dE 

dt 


jc^aV) 



drdTTr^pl^ioir). 


( 6 ) 


r 2 oo is the radius at which the mean density enclosed 
equals 200 times the cosmological mean at the redshift 
of formation of the halo: 

^= 200po[l + ^f(M)]®, (7) 

^^'200 

where po is the mean dark matter density at the present 
epoch (z = 0), and Zf(M) is the formation redshift of a 
halo of mass M. We may parameterize the halo density 
profile as: 


p{x) 


Ps 

a;“(I -I- x)l^ ’ 


( 8 ) 


The normalization constant A is chosen such that as = 
0.8, where as is the root mean square mass fluctuation 
in a sphere of radius 8 Mpc//i. The linear growth func¬ 
tion of matter overdensities is denoted by D{z). With 


where a; = r/rg is a dimensionless radius, pg and rg are 
constants for a halo. The well known Navarro-Frenk- 
White (NEW) |33j form is obtained when we set a = I 
and /3 = 2. 
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FIG. 1: We plot the dimensionless matter power spectrum for the standard cosmology (solid, black), as well as for the cases 
nis = 1.1, 1.2, and 1.5, for fcp = 100 h/Mpc (Panel a). An exponential cut-off is applied at the free streaming scale fcfs chosen 
to be 10® h/Mpc. Panel (b) shows the corresponding standard deviation cr of fluctuations normalized to ag = 0.8. 


(a) (b) 




FIG. 2: We plot the filling fraction /ah for different values of rus, for kp — 100 /i/Mpc (Panel a). The difference appears small 
at 2 = 0, but is substantial for large 2 . Panel (b) shows the redshift 2 * below which the halo contribution exceeds the free 
particle contribution. It is not very sensitive to the concentration parameter because of the exponential decrease of fm with 2 . 


We define the concentration parameter as: Note that we have defined r 2 oo and C 200 at the formation 

epoch. We may now express ps and in terms of M and 

^200 

C 200 — -■ 

rs 


( 9 ) 
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C 200 - Then the rate of energy release (Eq. is 
dEhalo {(Ta.v) 200 


dt 


rriy 


-M Pq[1 + Z({M)] /conc{c20o)) (10) 


where /conc(c 2 oo) is calculated for the density profile with 
concentration parameter C 200 as 


fconc — 


[/;^™dxa;2-«(l+x)-/3]" 


( 11 ) 


The lower limit e is required when the index a > 1.5. For 
the NEW profile with a = 1 and /3 = 2, we can set e = 0, 
and then Eq. is integrated analytically to yield 


fconc — 


^200 


1 — (1 + C200) 


-3 


3 [ln(l + C200) — C2oo(l + C200) 


2 • 


( 12 ) 


A halo of mass ~ IO^^Mq similar to the Milky Way is 
expected to have a concentration parameter C200 10 

[5^ . Earth-mass microhalos, on the other hand, are not 
expected to have large concentration parameters. [35] 
found concentration parameters C 200 ^ 3 for such very 
small halos. We simply assume a constant C200 = 5 inde¬ 
pendent of mass. 

To evaluate Eq. [T^ we need to determine the redshift 
of formation of the halo. Following the Press-Schechter 
formalism |36j , we assume that the probability of finding 
a halo of mass M at a redshift 2; is oc exp — [S^/2a^{M, z)], 
where 6c = 1.686 is the threshold for halo formation in 
linear theory. We can then estimate the averaged quan¬ 
tity: 


1 + Zf(M) 
1 + z 


3 4 f 

^ xi J a 


da; 


( 13 ) 
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where a;* = 6cl'j2a{z, M). Note that Eq. 
unity for large halo masses and large redshifts. 

The total energy due to WIMP annihilation per unit 
time and per unit volume may be obtained by integrating 


Eq. 10 over the halo distribution: 


of matter in bound halos: 

1 

h\\{z) = — - 

Po dM„i„ dM 

'ScD{0){l + z) 


dMM 


diV 


= erfc 


V2a{0,M) J ’ 

We also define the quantify ^(z) by 


(16) 


r°° dN 

/ dMM—{ 

■l + Zf(M)' 

] M ■ d-M 

[ 1 + z \ 


1 3 


)■ 

(17) 

which is larger than 1 at low redshifts but approaches 1 
for large z. Fig.j^a) shows the filling fraction for different 
choices of mg. We set kp = 100 /i/Mpc as our fiducial 
model parameter. The black curve shows the case rUg = 
Us = 0.96, i.e. the standard power law. Clearly, the 
filling factor can be many orders of magnitude larger at 
high redshifts if mg > Ug. 

It is also interesting to compute z*, the redshift below 
which the nonlinear halo contribution exceeds the free 
particle contribution. Fig. [^b) shows z* calculated for 
different values of the concentration parameter.lt is not 
very sensitive to the concentration parameter due to the 
exponential decrease of the mass function. For the stan¬ 
dard matter power spectrum, halos are only important at 
redshifts z ^ 50. However, for rUg = 1.5, halos contribute 
significantly even at z = 300. 

At high redshifts, we may make the “on the spot” ap¬ 
proximation when calculating the net energy input to 
the gas. Namely, we may safely ignore the propagation 
of high energy annihilation products from the redshift of 
emission to the redshift of absorption. The net energy 
absorbed per atom per unit time at a redshift z is given 
by: 


C(4 = 


/ dE 
nh{z) dtdE 

fWnV)rh /pcrit\ ^3 

ruy \ 


dE 

dtdV 



X 

1 - fm{z) + ^/conc(c20o)/fill(z)C(z) 

dM dt 


L 0 j 


(cTau) 200po 


fconc (C200)(l + zY 


r°° dN 

/ dMM^i 

■l-kZf(M)' 

M ■ dM ^ 

[ 1 + z J 


) (14) 


The comoving number density of halos is calculated 
from, for instance, the Press-Schecter mass function as 


diV 

dM 


Po 


1 dtJ 

udM 


exp 


' 2 cr 


1 2 


(15) 


We have simplified the mass function by setting po equal 
to the dark matter density, rather than the matter den¬ 
sity. Let us define the filling factor fmiz) as the fraction 


(18) 


where / is the fraction of energy absorbed by the gas, 
rib( 2 ) is the baryon number density, and fh is the mean 
nucleon mass, assuming 76% hydrogen and 24% helium. 
Hb, and are the baryon, dark matter, and total 
matter fractions at the present epoch. 


III. CMB CONSTRAINT 

In the previous section, we computed the rate of en¬ 
ergy release due to dark matter annihilation in halos. 
The fraction / of the released energy is absorbed by the 
gas, and a fraction »7ion(a:ion) of this energy goes into ion¬ 
ization, whereas a fraction ?7heat(a:ion) goes into heating. 
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FIG. 3: The excess optical depth due to dark matter annihilation inclnding the contribution from nonlinear haloes, for 
rria = 1.1,1.3,1.5 (Panel a). The magenta line shows At — 0.026 which is the maximum allowed excess optical depth above 
z = 6 (see text). Panel (b) shows the corresponding excess Compton y parameter due to dark matter annihilation as a function 
of z. 


We use the results obtained by m to estimate 77ion and 
Vhea,t- The ionization and temperature evolution follow 
the equations HU: 

- (1 + = ^J,[1 - Xioniz)]riioniz)^iz) 

dz 

- n{z)xl^iz)a(z) 

-il + z)H{z)^^ = -2T{z)Hiz)+^^^^^^az) 
, Xioniz)[T-^{z) -T{z)] 

+ W) ■ ^ ^ 

H « 0.07 eV“^ is the inverse of the average ionization 
energy per atom, neglecting double ionization of helium. 
In the above equations, a is the case-B recombination co¬ 
efficient, Try is the CMB temperature, fcb is Boltzmann’s 
constant, and is the Compton cooling time scale « 
1.44 Myr [30/(1-|-x)]^. The last term in the temperature 
evolution equation accounts for the transfer of energy be¬ 
tween free electrons and the CMB by Compton scatter¬ 
ing |38H40) . In the temperature coupling term, we assume 
a^ion 1 and ignore the helium number fraction. In prac¬ 
tice, we compute xion and Tgas using a modified version 
of the publicly available recfast program [3^ [40] . 

CMB photons scatter off free electrons that are present 
due to partial ionization of the gas. Thomson scat¬ 
tering of the CMB causes damping of the temperature 
anisotropy TT power spectrum, as well as a boost in the 
large angle EE polarization power spectrum |16j . The 
scattering is quantified by means of the optical depth de¬ 


fined as the scattering cross section times the free electron 
density integrated along the line of sight: 

t{zi,z) = / dtcarneiz) (20) 


where ctt is the Thomson cross section. We also calculate 
the ‘excess’ contribution as 

CCTt k crit 

^ Hioom VfW? 

X f dz{l + Ax{z) (21) 

Jzi 

with i7ioo = 100 km/s/Mpc and Ax = Xion — a^stdj 
where Xstd denotes the standard recombination history 
(i.e. without dark matter annihilation, and rus = ris). 
We have ignored dark energy, and therefore, the above 
equation holds true for zi^^ d. 

We see that even small changes in xion at high redshifts 
can boost the total optical depth due to the ^/1 -I- z term. 
One may also hope to constrain dark matter annihilation 
by measuring the spectral distortion of the CMB, quan¬ 
tified by the Compton y parameter [mill]: 


y = 


y - 2 /std « 


X 


fcb [T{z)-Try{z)] 

J rUeC^ 

CCTt k crit /Ifl Oblt^ fcb 
Hloo m meC^ 


[ dz{l + zy/^ [x^tdAT + Ax{T,td - Try)], (22) 

J Z-\ 
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FIG. 4: The top panels (a) and (b) show the ionization and temperature history of the Universe. The black curve is plotted 
for the standard ACDM cosmology, i.e. ignoring dark matter annihilation. The red curve {rus = Us) is for a standard power 
spectrum, but includes the effect of dark matter annihilation. The green, blue, and magenta curves are plotted for the modified 
power law of Eq. The bottom panels (c) and (d) show the temperature and polarization power spectra for the four models. 


where AT(z) = T{z) — T'std('Z); and as before, T^td rep¬ 
resents the gas temperature in the standard ACDM sce¬ 
nario. 

Fig. I^a) shows the excess optical depth due to dark 
matter annihilation with rUg = 1.1, 1.3, and 1.5, with 
dark matter mass = 100 GeV, and concentration pa¬ 
rameter C 200 = 5. Assuming the total optical depth mea¬ 
sured by Planck r « 0.066 m and full ionization up to 
z = 6, we find that the excess optical depth At < 0.026, 


which excludes large values of rUg. The Compton y pa¬ 
rameter is less constraining because it is weighted to¬ 
wards large z when the gas temperature is close to the 
CMB temperature due to efficient Compton scattering. 
The planned pixie mission can constrain |i/| < 2 x 10“® 
|42) . and may exclude very large values of rUg at the rel¬ 
evant length-scales. 

Fig. i shows the evolution of the ionization fraction 
(Panel a) and the gas temperature (Panel b), for dif- 
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ferent values of rUs- We assume f{as,v)jm^ = 1/100 
pbxc/GeV for the figure. The black curve is plotted for 
the standard ACDM, i.e. ignoring dark matter annihi¬ 
lation. The red curve is plotted for the standard power 
law, but with accounting for dark matter annihilation. 
The green, blue, and magenta curves are the results for 
TOg = 1.1, 1.3, and 1.5, respectively. The effect on the 
CMB TT and EE power spectra is shown in Panels (c) 
and (d). For large mg, significant damping is caused on 
the TT power spectrum, as clearly seen in Panel (c). 
Also plotted in (c) are the TT power measurement from 
Planck. It might appear that large values of Wg are al¬ 
ready excluded at high significance. However, the am¬ 
plitude of the CMB power spectrum is determined by 
^ Agexp—2r. While the optical depth r is increased 
by ionization by dark matter annihilation, the effect is 
almost fully degenerate with the amplitude of the pri¬ 
mordial curvature power spectrum Ag, except on very 
large scales that were outside the horizon at the time of 
particle annihilation. Therefore, the TT power spectrum 
alone cannot be used to place constraints on dark mat¬ 
ter annihilation, but the degeneracy is broken by using 
information of the CMB polarization, because Thomson 
scattering causes a boost in the large angle polarization 
power spectrum. A second technique to break the degen¬ 
eracy is through the measurement of gravitational leasing 
of the CMB by large scale structure. The Planck exper¬ 
iment has recently measured the lensing potential at the 
40(7 level [43]. Measurement of the gravitational lens¬ 
ing of the CMB places constraints on the combination 
= 0.591 ± 0.021. Since other constraints exist 
for Hm, the measurement of gravitational lensing places 
a bound on erg, and hence on Ag. Gravitational lensing 
thus breaks the degeneracy between Ag and r. The recent 
results from Planck give us log(10^°Ag) = 3.064 ± 0.023, 
and T — 0.066±0.012, which we can now be used to place 
bounds on the root mean square mass fluctuation (imax, 
and hence on the power law index rUg. 

We consider a number of models with different mg 
and calculate (Tmax- We also compute the correspond¬ 
ing ionization history and the CMB power spectra us¬ 
ing a modified version of the game software, assum¬ 
ing pann = /(^au)/??!^ = (1/100) pbxc/GeV. Foi' each 
model, we fit the theoretical power spectra to the ob¬ 
served TT power spectrum, by modifying the quantity 
AgC”^"^. By means of Montecarlo simulations, we obtain 
a bound on the combination Ase”^"^ = 1.872 ± 0.101. 
The 2(7 upper bound on AgC”^’’ results in a correspond¬ 
ing bound on the standard deviation of mass fluctuations: 
CTmax < 100 at the 2a level. This finally translates to a 
bound on the power law index: rUg < 1.43(1.63) for kp = 
100 (1000) h/Mpe. 

IV. CONCLUSIONS 

We have proposed a new probe of the matter power 
spectrum on very small scales, through the effect of dark 


matter annihilation on the thermal evolution and on the 
ionization history of the inter-galactic gas. We have 
considered a simple modification to the standard power 
law, for the primordial power spectrum: Pprim(fc) 
(A:/fcp)'"% where rUg > rig. Such a power law is accept¬ 
able provided kp is large enough, e.g. kp > 10 h/Mpe. 
The form of the root mean square mass fluctuation cr has 
been calculated as a function of the power law index mg. 
The maximum value a^a.^ = o’(Mmin) varies significantly 
with mg. One may expect a 1 — cr fluctuation to enter the 
non-linear regime when 1 ± z « (7max/<5cP(0)- We have 
then computed the filling fraction (the fraction of dark 
matter bound in nonlinear halos) as a function of red- 
shift. For large mg, there are many orders of magnitude 
more halos at z > 100. For the standard power spectrum 
power law, halos are only important for z ^ 50. On the 
other hand, when mg ^ 1.5, halos provide the dominant 
contribution to total dark matter annihilation rate even 
at z = 300. 

We derived an explicit expression for the energy in¬ 
jected per unit gas atom per unit time at a redshift z. A 
large contribution from dark matter halos at z ~ 300 can 
significantly alter the spectrum of CMB anisotropies. We 
computed the CMB power spectra using the came code, 
for a fiducial dark matter annihilation cross-section of 
Pann = /(CTaw)/m^ = (1/100) pbxc/GeV. We used the 
Planck (2015) TT power spectrum data to test theoret¬ 
ical models. The current data already excludes a root 
mean square fluctuation (7niax = cr(Mi„in) ^ 100 at the 
2(7 level. The bound on (7max may be expressed as a 
constraint on the power law index on small scales: We 
exclude mg > 1.43(1.63) for kp = 100 (1000) h/Mpe. 
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